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Introduction: The Lya absorption line forest, seen in quasar spectra, is generally interpreted as 
being due to cosmologically distributed “clouds” of primordial gas. Analyses of the observations {e.g. Rauch 
ei al. 1992) reveal that the number distribution can be described by power laws: dAf/dz oc (1 + zf and 
dAf/dNffi oc , where M hi is the HI column density. The typical values for power law indices range 
between 2 < 7 < 2.6 and 1.7 £ /? < 1.9. One model postulates that the Lya clouds are optically thin 
entities, photoionized by the background UV flux, J„ a (1 + z) 7 , and confined by an adiabatically evolving 
intercloud medium (ICM): P{z) oc (1 + z) 5 . Analytic studies of this model suggest that the ensuing Lya line 
statistics can account for the observations (in particular, the dM/dz and the dM/dNm distributions) if the 
cloud mass spectrum is a power law dM JdM oc M~ s , 6 ss 1.9, and j ft s 4 ( c.f Ikeuchi k Ostriker 1986, 10). 
One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for 
the clouds at all epochs, the validity of which is questionable. 

Simulations: We investigate the pressure-confined model using a 1-D spherically symmetric 
hydrodynamical code to simulate cloud evolution over the epoch 1.8 < z < 6 This enables us to re- 
lax many of the assumptions incorporated in the analytic studies. We only consider clouds with masses 
2.75 < log (M/ M©) < 9.25. Clouds with larger masses are Jeans unstable and tend to collapse; they are not 
“pressure-confined” . The lower mass limit is chosen to ensure all clouds that can produce absorption features 
with equivalent widths W r > 0.2A over the epoch 1.8 < z < 3.6 are considered. We calculate Voigt profiles 
from the model clouds at varying impact parameters and use the resulting absorption features to determine 
line-profile parameters {e.g. Nhi, the rest-frame equivalent width W r , and the Doppler parameter b). We 
construct synthetic samples of lines subject to similar selection criteria used to define a sample of ~ 400 
observed lines (i.e. W r > 0.2A). In this sense, the present study is unique: past theoretical studies have 
relied on column density thresholds. Nhi , however, is a fitted parameter and is sometimes subject to large 
uncertainties. Line widths have been traditionally used to define observational samples as they are directly 
measurable. We compare the synthetic dM/dW r , dN/dNni and dAf/dz distributions against the observed 
trends. The synthetic results are normalized by demanding that the number of clouds with W r > 0.2A be 
the same as observed. Apart from studying the properties of clouds immersed in j = 4 UV background, we 
also consider the possibility that the UV background is constant at the epoch of interest. For more details, 
see Williger k Babul (1992; WB). 

Results: The redshift-integrated dAf /dW r distributions for j = 0, 4 (5 = 1.90) are shown in Figure 
1. Clearly, there is a shortage of model lines with W r > 0.3A, in comparison with the observations. The 
result is partly due to the fact that (6) ss 22kms -1 for simulated pressure-confined clouds (with a dispersion 
of approx, lkms -1 ) — i.e. b is due to thermal velocities — while (6) ss 35kms -1 for the observed lines 
(with a dispersion of approx. 15kms -1 ). The lower 6 for the simulated clouds implies that these clouds 
require larger column densities than those associated with the observed lines in order to account for the 
observed dAf /dW r distribution. The imposition of the upper mass limit for the pressure-confined clouds, 
however, restricts the column density range of the simulated clouds. 

The dM/dNni distributions for the simulated and the observed clouds in the redshift range 2.6 < 
z < 2.8 (Figure 2) shows that although the two are consistent with each other for log(Vj//) 15.2, there are 
fewer model clouds with Nhi larger than observed. Note that the shortage is more severe for the j = 4 case 
as the UV background at z « 2.7 is more intense than in the j = 0 case and therefore the clouds are more 
highly ionized. We also show the “analytic” dAf /dN hi, computed using the simple homologously-expanding 
pressure-confined cloud model of 10 and the prescription of WB {e.g. we explicitly take into account the upper 
mass limit). This distribution cuts off at even a lower Nhi than the synthetic distribution; the simulations 
reveal that the more massive clouds tend not to expand homogolously and, consequently, develop non-uniform 
density profiles. Self-gravity also retards their expansion. 

In Figure 3, we show the observed and the synthetic dAf /dz distributions for clouds with \V r > 0.2A; 
the dAf /dz distribution for a subset of clouds with W r > 0.3A is shown in Figure 4. For W r > 0.2A sample, 
the synthetic results match the observations well. Neither the normalization (reflecting the discrepancy in 
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the dffjdWr distributions) nor the slope of the synthetic dAf/dz (the number of pressure-confined clouds 
decreases much too rapidly at low redshifts) matches the observations for the subset with W r > 0.3A. The 
rapid decline in clouds with large equivalent widths towards lower redshifts is, once again, due to the upper 
mass limit for the clouds. 

Conclusions: As a result of an upper limit to the mass of pressure-confined clouds and the 
“thermal” Doppler parameters, the equivalent width and column density distributions for spherical, pressure- 
confined clouds with a power-law mass spectrum fail to account for the observations. The upper mass limit is 
determined by nature of the evolution of both the ICM pressure and the UV background; in order to achieve 
an acceptable fit to the data, we require the ICM pressure to decrease as P(z) oc (1 + z) p , p > 5.8 (WB). 
This is difficult to justify physically. Aspherical clouds, on the other hand, may allow the pressure-confined 
model to remain viable. 
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Fig. 1. dtf/dW r of the observed clouds vs. the z — integrated 
distribution of the model clouds. The j = 0, 4 cases for 6 = 
1.90 are shown by solid, dashed lines respectively. The error bars 
indicate bin widths and lcr errors for the observations. 


Fig. 2. Column density distribution of the observed clouds vs. 
models in the redshift range 2.6 < z .< 2.8. The various curves 
correspond to various analytic and synthetic results: j = 0/j = 
4 synthetic (solid/dashed); j = 0/j = 4 analytic (dotted/dot- 
dashed). At high N Hl , the synthetic curves steepen before an 
abrupt cutoff. The cutoff is even lower for the analytic models. 




Figs. 3-4. The observed vs. model redshift distribution of Lyo clouds with W r > 0.2A (Fig. 3) and 
W r > 0.3A (Fig. 4). The various curves correspond to various analytic and synthetic results: j = 0/j = 4 
synthetic (solid/dashed); j = 0/j = 4 analytic (dotted/dot-dashed). Although the synthetic curves fit the 
data for W r > 0,2 A, there are significant discrepancies in the slope and the normalization for W r > 0.3 A. 
The analytic fits are worse, except for j — 4, W r > 0.2A. 
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